A CENTRAL LIMIT THEOREM FOR THE SPECTRUM 
OF THE MODULAR DOMAIN 
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ZEEV RUDNICK 

The statistics of the high-lying eigenvalues of the Laplacian on a Rie- 
mannian manifold have been intensively studied in the past few years 
by physicists working on "quantum chaos" . A number of fundamental 
insights have emerged from these studies, though to date these have 
yet to be set on rigorous footing. In the case the manifold at hand is of 
arithmetic origin, these studies are related to some profound number 
theoretical problems and as such may be more amenable to investiga- 
tion. In this note I make use of the arithmetic structure of the modular 
domain to establish Gaussian fluctuations in its spectrum for certain 
smooth counting functions. 
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w{Ej < x\ r^ — — X, as x ^ cxD 



1. Background 

To set the stage, I start with describing some of what is currently 

believed to hold for the statistics of the eigenvalues. Given a compact 

Riemannian surface M, Weyl's law for the eigenvalues Ej of the Lapla- 

cian says that the number of eigenvalues below x grows linearly with 

x: 

area(M) 

47r 

Let n(E, C) be the number of levels in a window around E for which 
the leading term in Weyl's law predicts C levels: 

27r 2tt 

^ ' ^^ area(M) ^ area(M) ^ 

and more generally for a test function / define 

(1.1) n,(E,/:):=V/("^^"^^) (^^-^^^ 



^^^ Atx C 



J 



which counts the levels lying in a "soft" window of length — HttkC 

JO o area(Af) 

about E. In the above C = C{E) depends on the location E. In what 
follows we will usually write n{C) for n{E,C), the dependence on E 
implicitly understood. 

To study the statistical behaviour of n{C) we need to consider E as 
random, drawn from a certain distribution on the line. We denote by 
(■) this kind of energy averaging, e.g. {F) = ^ f^ F{E')dE'. Weyl's 
law leads us to expect that the mean value of n{C) is C and likewise 
that of nf[C) is £ ■ / f{x)dx. 

1.1. Number variance. The variance of nf{E,C) from its expected 
value is: 

E}{E,C) = {\nj{C)-MC))\') 
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It is customary to express the number variance by means of an integral 
kernel Ke{t), called the "form factor", so that as E —>■ oo 

Sj(E, C)^C- f{ufKE{-^)du = / {Cf{CT)fKE{T)dT . 

where f{y) = j_ f{x)e~^'^^^ydx is the Fourier transform of /. 

For "generic" surfaces, Berry j31 Ej argued that as i^ -^ oo, the 
behaviour of ^"jiE, C) for C in the range^ 



(1.2) l^^C^^ Cmax = VE 

is universal, depending only on the coarse dynamical nature of the 
geodesic flow on the surface, and follows that of one of a small number 
of random matrix ensembles: If the flow is integrable (as in the case of 
a flat torus) then T,'j{E,C) ~ C ■ /^ f{x)^dx for £ — > oo, as in the 
Poisson model of uncorrelated levels. If the flow is chaotic (as in the 
case of negative curvature) then the behaviour is as in the Gaussian 
Orthogonal Ensemble (GOE): For the sharp window (/ = l[-i/2,i/2])5 
this is given by T?{E,C) ~ ^log£ for £ — ;► oo. For sufficiently 
smooth /, in the GOE we have T?AE,C) ~ 2 j'^^f{uf'\u\du, that is 
the variance of sufficiently smooth statistics tends to a finite value as 
£ — > oo. The form factors for the random models are K^°'^^{t) = 1, 
and 



GOE.,^^^2|r|-|r|log(l + 2|r|), |r| < 1 
2-|r|logi±^, |r|>l 



K---{T 



It is to be emphasized that the above behaviour is only valid in the 
universal regime 1 ^^ £ ^^ ^/E; for C yy y/E the integrable case 
is fairly well understood (at a rigorous level), see the survey f5^: The 
variance grows as a/E (a classical result [T2j in the case of the standard 
flat torus). In the chaotic case it is believed [3111 that generically, the 
number variance continues to be small as in the universal regime. 

1.2. Fluctuations. Our main interest here is in the value distribution 
of the normalized linear statistic 

nf{E, C)-{nf{C)) 



as E varies. In all the statistical models (Poisson and GOE/GUE), it 
is a standard Gaussian j23IIlllini- 



^The symbol f{x) ^^ g{x) means that f{x)/g{x) 
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In the integrable case, when L y>- yE, the distribution is known 
([T7j. [5]), and is definitely not Gaussian. Inside the universal regime 
(|1.2|) . the distribution is believed to be Gaussian in both the integrable 
[0] and chaotic PJ IHl ES] cases. In the special case of the standard fiat 
torus, this has been proved in a small part of the universal regime near 

/:e 1181. 



1.3. The modular domain. We start with the upper half-plane EI = 

{x + y/—ly : y > 0} equipped with the hyperbolic metric ds"^ = 
y~'^{dx^ + dy"^), which has constant curvature equal —1. The Laplace- 
Beltrami operator for this metric is given by A = y'^{-§^ + ■§-^)- The 
orientation-preserving isometrics of the metric ds^ are the linear frac- 
tional transformations P^L(2,M) = 5L(2, M)/{±J}. 

The modular domain is the Riemann surface obtained by identifying 
points in the upper half plane which differ by a linear fractional trans- 
formation with integer coefficients, that is by elements of the modular 
group 

r := PSL{2,Z) = SL{2,Z)/{±I} . 

The resulting surface H/F is non-compact and has cone points, but has 
finite hyperbolic area: area(]HI/r) = tt/3. 

The spectrum of the Laplacian on the modular domain has a con- 
tinuous component. Nonetheless, Selberg j21] showed that a version of 
Weyl's law holds for the discrete spectrum: If we write the eigenvalues 
of -A on L2(H/r) in the form Ej = 1/4 + rj, then: 

(1.3) #{r, <T}= ^^"<W^) t^ - ^TlogT + c^T + 0( ^ ^ 



477 TT log T 



where Ci = — — 



In the case of the modular domain, deviations from generic statis- 
tics were discovered JTUIIZIITT]. Although the geodesic fiow is chaotic, 
the local statistics of the spectrum seem Poissonian, and Bogomolny, 
Leyvraz and Schmit jS] argued that the behaviour of the form factor is 
given by 

{ exp(c2VEr) _2 ^ , , , \ogVE 

for some constants Ci, C2 > 0, that is to say, in the universal regime we 
have 

y2fj, ,. J2c,^J-f{ureMc2^)du, ^^^C^^VE 
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The only rigorous results known concern the the closely related case 
where the modular group is replaced by quaternion groups: In the 
1970's Selberg [13 Chapter 2.18] gave a lower hound for the variance 
S^(£', E) of the sharp counting function n{E) of the form E^(£', E) 3> 
\^/{\ogEy. Luo and Sarnak |20j gave lower bounds for the averaged 
number variance of the sharp counting function n[E,C): 

-/ J:\E,C')dC' :^Ve/ {log Ef, \fE / \og E ^^ C^^\fE 
>-' Jo 

in the case of arithmetic (co-compact) groups. In the case of the hard 
window / = l[-i/2,i/2] no upper bounds are currently available. 



2. Formulation of results 

2.1. Definition of the smooth counting functions. Let / be an 

even test function, whose Fourier transform / G C^(]R) is smooth and 
compactly supported, and normalized by requiring that 



f{x)dx = 1, 

and that 

sup{|e|:/(O^0} = l. 

The eigenvalues of the Laplacian are parametrized by Ej = 1/4 + r|. 
We define smooth counting functions by 

(2.1) Nf,L{T) = Y, IWr, - r)) + /(L(-r, - r)) 

i>o 

This in essence carries the same information as (jl.l|) . The relation 
between the expected number of levels C and the inverse width L of 
the momentum window is 



area(M/r)V^ _ VE 
The leading order behaviour of A^/,l is given by 



NfAr) ■■= ^ I {f{L{r - r)) + /(L(-r - r))}M(r)dr 



where 



area(H/r) , , , F', , F' 1 

M{r) := ^--^^rtanh(7rr) - —(1 + ir) - — (- + ir) . 

^ J. J. ^ 
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(In keeping with tradition, I use the symbol F for both the modular 
group and the Gamma function). The term Nf^^ij-) is asymptotic to 
£: 

area(e/r) ^ r ^.logr Ir 



-—^ , ^area(e/r) r ,, ^, t ^, 
^fAr) ~ 2 A_^^ y /(a;)da;- + 0{ 



L ' 6L 



2.2. The results. We will see that A^/^i—A^/^i has mean zero and show 
that the variance oi Nf^^, when lim sup ttL/ log T < 1, is asymptotic to 



(2.2) al := ^ 1^^ fiufe-^-du 

where k = ^ Up^2i'^ + ^^^^) = 1-328 .... Thus when the expected 
number of levels C satisfies 



area(H/r)-^-^ < C ^^ T = Ve . 
logE 

the form factor Ke{t) is given by 

expc2\/^r 
^bM = ci — 



with Ci = Gk/tt, C2 = 7r/6. 

Our main result is that the fluctuations of A^/,l are Gaussian: 

Theorem 2.1. Assume that L —>■ oo as T ^ oo but L = o(logT). 
Then the limiting value distribution of {Nf^^ — Nf^^j/ai is a standard 
Gaussian, that is 



1 ^.„.r^ r- irr ^rr., . NfA^) - Nf^r) ^ ^, _ f .-nV2^ 



lim - measjr G \T, 2T] : ^'^ ' ^^^^^^ < x| = / e ^ 

The reason that we need to assume that L = o(logT) is that we prove 
Theorem 12.11 by the method of moments, and in computing the i^-th 
moment we find Gaussian moments for L < ck log T, where c;^ — > as 
K grows. 

Plan of the paper: In sections 0] and 0] we give some results on the 
hyperbolic conjugacy classes of the modular group. In section El we use 
Selberg's trace formula to express Nf^iij) — Nf^iij) as a sum Sf^iij) 
over hyperbolic conjugacy classes plus a negligible term. The variance 
of Sf^L is computed in section IHl and the higher moments in section [7| 
We prove Theorem 12.11 in section |H1 
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3. The modular group 

3.1. Conjugacy classes. To analyze iV/,L we use the Selberg trace 
formula, which for a discrete co-finite subgroup F C PSL{2, M) relates 
a sum over the spectrum of the Laplacian on L^(EI/r) with a sum over 
the conjugacy classes of the group F. We review some background 
material on these classes for the modular group PSL{2, Z). 

The conjugacy classes are divided into the class which consists of 
the identity element, hyperbolic, elliptic and parabolic classes. The 
hyperbolic conjugacy classes in F are represented by matrices P which 
are diagonalizable over the reals and are conjugate to a matrix of the 

form I \-i I with A > 1. The norm of P is defined as A/'(P) = A^. 



A- 

The norm is therefore related to the trace of the corresponding group 
element by 

Af{py/^+M{P)-^/^ = |tr(P)| 

We can write as each such P as, P = Pq where Pq is primitive and 
/c > 1, where an element of F is primitive if it cannot be written as an 
essential power of another element. As is well known, primitive hyper- 
bolic conjugacy classes correspond to closed geodesies on the Riemann 
surface H/F. 

In the case of the modular group F = P5'L(2,Z), the traces are 
integers. If P is a hyperbolic class with trace | tr(P)| = n, n > 2 then 
its norm is 

(3,1) Min) = ( " + ^ f 

For the modular group, primitive hyperbolic conjugacy classes are 
parametrized by indefinite binary quadratic forms as follows (cf. j23j): 
Take a binary quadratic form Qa,b,c{x,y) := ax"^ + bxy + cy"^, with 
a,b,c E Z. The discriminant of Qa,b,c is d := 6^ — 4ac. The form Qa,b,c 
is indefinite iff rf > 0. We assume that d is not a perfect square. We say 
that Qa,b,c is primitive if gcd(a, 6, c) = 1. Two binary quadratic forms 
Q, Q' are equivalent if Q'{x, y) = Q{ax + by, ex + dy) for an element 

7=1 J I of SL(2,Z); since the forms are quadratic, they are also 

equivalent under —7 and hence equivalence is over PSL{2,Z). Let 
h{d) be the number of equivalence classes of primitive binary quadratic 
forms of discriminant d. 

The automorphs of Qa,b,c are all of the form 

±i-[t,u) ±1 ^^ Ut + bu) 
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where (t, u) solve the Pelhan equation 

(3.2) t^ -du^ = 4 

If n 7^ then these are hyperbohc elements of SL{2, Z) with norm 
^f{P) = {t + u\fd)'^ and trace t. 

Let e^ = |(td + Udvd) (t^, Ud > 0) be the fundamental solution of 
(j3.2j) . Then the matrix P{td,Ud) is a primitive hyperbolic matrix Pq 
of trace tr(Po) = td and norm Af{Po) = e^. It turns out that in this 
way we get a bijection between equivalence classes of primitive binary 
quadratic forms and conjugacy classes of primitive hyperbolic matrices 
in PSL{2,Z). Thus the number of primitive hyperbolic conjugacy 
classes of norm e^ is precisely the class number h{d). 

3.2. The amplitude P{n). We define, for n > 2, 

.^ ^^ ^f ^ _^ \- logAr(Po) 

l^-^^ PW-^ 2^ Ar(P)V2_Ar(P)-V2 

tr{P)=n ^ ' ^ ' 

the sum over all conjugacy classes {P} in PSL{2, Z) with | tr(P)| = n, 
equivalently with norm A/'(n) given by ()3.1|) . These quantities turn 
out to be crucial in our analysis. The factor 1/2 in the definition is 
inserted among other reasons to give numbers with mean value 1 (as 
can be seen from the Prime Geodesic Theorem): 

^ /?(n) ~ A^, as A^ ^ oo . 

n<N 

As representatives of the conjugacy classes of matrices with trace 
ra > 2 we can take the matrices Pq = P{td,Ud)^ = P{n,u) where d 
runs over all discriminants, k > 1 and n^ — du"^ = 4, n > 2, u > 1. 
Thus we see that 



(3.4) (3{n)= Y. 



h{d) log ed 



d,u>l v^*^ 



-du^ 



Dirichlet's class number formula allows us to use ()3.4j] to express l3{n) 
in terms of Dirichlet L-functions: For a discriminant d on associates 
the quadratic character Xd given by Xdip) = ( - ) for P an odd prime, 

Xd(2) = 1 if rf = 1 mods, Xd(2) = -1 if rf = 5 mod 8 and Xd(-l) = 1- 
The associated L-function is L{s,Xd) = J2n>iXd{n)n~'^, 9^e(s) > 1. 
Dirichlet's class number formula is 

h{d) log ed = ^/dL{l,Xd) ■ 
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Inserting the class number formula into ()H.4|1 we find that 

(3.5) /?(n)= J2 l^il^Xd). 

d,u>l:dv?='n?—A 

As a consequence, one can get an upper bound of /3(n) = 0((logn)^) 
by using L(l,Xd) <C \ogd. What is more useful to us is that, in the 
mean square, /3(n) is constant: 

Lemma 3.1 (M. Peter [21j). 

(3.6) ^ ^{nf ^ kN, N ^oo 

n<N 

where k, is given by the product over primes 
, „, 1015 TT, n^-2)o3 + l 

This (complicated) expression was derived heuristically by Bogo- 
molny, Leyvraz and Schmit [8^ and proven by Manfred Peter J2Tj, who 
uses the expression (13. 5j) for /3(n) in terms of L(l, x) and methods re- 
lated to work on moments of class numbers P^. For an extension to 
the case of congruence groups, see J19j. 

4. The length spectrum 
We will need to study alternating sums of the form 

K 

5^±logAr(n,) 
i=i 
The first question is to when these alternating sums vanish. 
We say that a relation 

K 

(4.1) ^r7,logAr(n,) = 0, r/, = ±1 

i=i 
is non- degenerate if no sub-sum vanishes, that is if there is no proper 
subset 5* C {1, . . . , K^ for which ^.g^^/j logA/'(raj) = 0. The existence 
of non-degenerate relations ()4.H) forces severe constraints. To explain 
these, recall that M{n) is a unit in the real quadratic field Q(a/'^^ — 4). 
We claim that such such relations can occur only if all these units lie 
in the same quadratic field. 

Lemma 4.1. Let ^ • ±logA/'(nj) = be a non-degenerate relation. 
Then all the norms N'ini) lie in the same quadratic field, that is for 
some common d we have n^ — 4 = dff' for all i. 
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Proof. We can write each norm as a power of the fundamental unit of 
the quadratic field in which it lies. Thus it will suffice to show that if 
Fi, . . . , Fk be distinct real quadratic fields, then the fundamental units 
Si of Fi are multiplicatively independent. 

Let E = FiV ■ ■ -V Fk be the compositum of the fields F^. This is a 
Galois extension of the rationals with Galois group G = Gal{E/Q) an 
elementary Abelian 2-group (Z/2Z)*, for some s < K. If we denote by 
Ue the unit group of E, then G acts on Ue and hence we get a linear 
representation on the vector space Q,®Ue- 

We claim that the e^ are eigenvectors of G, that is cr{ei) = ef for 
all a G G, where Xi '■ G -^ {±1} are distinct characters. This forces 
them to be multiplicatively independent. 

Indeed, since we have an Abelian extension, all subfields are Galois 
and in particular Fi are preserved by G. Since the unit group is also 
preserved this means that under the action of any element a G G, ej 
is taken to a unit of Fi which is necessarily e^^. That is we have a 
character Xi of G with o'(ej) = e^ . The characters Xi are distinct 
since the kernel of Xi is precisely Gal{E/Fi). D 



We next get a lower bound for ^ . ^ ±logA/'(nj) in the case it is 
non-zero. 

Lemma 4.2. i) If m ^ n then 

I logJ\f{m) - logA/'(n)| > 



mm(m, n) 



K 



a) Suppose X]7=i ='=logA/'(nj) is nonzero. Then 

1 



J]±logAr(n,)|>- 

(nf=iAr(n,; 



2^-1-1/2 



Proof, i) Indeed, since logA/'(n) = 21ogn + 0(l/n^), if m 7^ n, say 
m > n, then 

771 1 

logA/'(m) -logA/'(n) = 21og hO' 



n n"^' 



Since log ^ > log ^^ > i, we find 

1 

log Af{m) - \ogM{n) > - 



n mm m,n 
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ii) Let Xj = "'^\^ so that Af{nj) = X], and set a = H^i >^f^- If 
|a- 1| < 1/2 then 

K 

I ^±logA/'(nj)| = 2|loga| > |a - 1| 
i=i 

So it suffice to give a lower bound for |q; — 1|, assuming a ^ 1. 
This follows from Liouville's theorem on Diophantine approximation 
of algebraic numbers by rationals; we give an explicit proof as follows: 
Let E = Q(y^, . . . , y/uK), which is a Galois extension of the rationals 
with Galois group G = Gal{E/Q) which is an elementary abelian 2- 
group of order 2*, for some s < K. Moreover a G ii^ is an algebraic 
integer, and hence the norm A/'s/(Q(a; — 1) is a nonzero rational integer, 
hence has absolute value at least 1. Thus 

\AfE/Q{a - 1)\ = \a - 1\ JJ |a'^-l|>l 
Since AJ = A^^ for all a E G, we have 

K 

la'^-ll^HX. + l 
i=i 



Thus 



|a-l| > ,_^ — ^-> 



(n;ii A, + 1)1^1-1 (n;=iArK))(^"-^)/^ 

D 

5. An EXPANSION FOR Nf^L 

5.1. The Selberg Trace Formula. We will transform Nf^i by using 
the Selberg trace formula j21|: Let g G C^(M) be a an even, smooth 
and compactly supported function, and let 



h{r) = I g{u)e"'''du 
so that 



giu) = ^ r /i(r)e-™dr . 
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The Selberg Trace Formula for a discrete co-compact sub-group F C 
PSL{2, M) with no elhptic elements is the identity [23 

{P} hyperbolic ^ ^ ^ ^ 

where the sum is over all hyperbolic conjugacy classes of F. 

In the case of the modular group, the hyperbolic terms can be written 
as 

(5.2) 2-J2Pin)9{iogAf{n)) 

n>2 

where the amplitude [3{n) is given by ()3.3|) . 

For groups with elliptic elements, there is an extra contribution to 
the RHS of ()5.ip which is a sum over the finitely many conjugacy classes 
of elements E of finite order m > 2: 

"^ 1 /"oo —2nkr/m 

(5.3) yy . , , , , / h{r)^ —dr 

For discrete groups whose fundamental domain is non-compact but 
of finite volume, that is with cusps, there are extra terms coming from 
the contribution of the continuous spectrum and parabolic elements. 
For F = PSL{2, Z), these terms are given explicitly by [T6] : 

(5.4) g{0)\og--—J /i(r)(-(l + ir) + -(- + ir)Jdr 

+ 2V-^^(21ogn) 
^-^ n 

n=l 

where A(n) is the von Mangold function. 

5.2. Transforming A^/,l- We now apply the trace formula to derive 
an alternative expression for A^/,l- Taking 

h{r) = f{L{r-T)) + f{L{-r-r)) 



so that 



we find that 



«(■")- 2^^'(2^)(^""" + ^'"') 



(5,5) N,,dT) = A'/,i(r) + S;,i(T) + S 

where: 
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• The term A^/,l is given by the contribution of the identity class to 
()5.H) and part of the paraboHc terms in ()5.4|) : 

1 



NfAr) = ^ / {f{L{r - r)) + f{L{-r - r))}M(r)dr 



where 

,^, , area(H/r) , , , T'^ . , T' 1 . , 

M{r) = ^-^^r tanh(7rr) - —(1 + ir) - — (- + ir) 

By Stirhng's formula, we have 

1 f°° ^, ,, T ^,logr. 



iV/,L(r) = - / f{x)dx- + Oi ^ 



— oo 



• The term S'/i(r) is the contribution of the hyperbolic classes (j5.2p : 

(5.6) Sf,L{r) = i- y /5(n)/(i^^^^) (e--i°sArW + g-iogA^W) 

ttL ^—^ 2tcL 

n>2 

The sum ()5.6|) contains only terms with logA/'(n) < 2'kL, that is n < 

As we will see below, it is the term Sf^ilr) which is responsible for 
the fluctuations of Nj^i^r), and its variance is asymptotic to cr|_. As we 

can see from the formula ()5.6p . since / has compact support we have 
Sf^Lir) = for L < 1. 

• S is the contribution of the elliptic classes ()5.3|) and the remaining 
part of the parabolic contribution ()5.4j) , namely 

(5.7) J-/(0)logl + i-f ^/(!2i!i)2cos(2rlog„). 

n=l 

£ is easily seen to be negligible, that is £^ = oi^aL). Indeed, the con- 
tribution of the elliptic elements is easily seen to be 0(e^*^*^^ /L). 
As for ()5.7|) . this is bounded as L ^ oo by (say) Mertens' theorem. 
Moreover the mean value of ()5.7jl clearly vanishes as T ^ cxd. 

We thus see that the difference between the centered counting func- 
tion Nf^ii'^) — Nf^L and the sum ^/^^(r) over hyperbolic conjugacy 
classes is negligible relative to the standard deviation ai of Sj^l^t), 
and thus for our purposes we need only investigate the statistics of 

6. The mean and variance of Sf^L 

6.1. The averaging procedure. We define an averaging procedure 
by taking a non- negative weight function w > 0, which is smooth and 
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compactly supported in (0, oo), with f^ w{x)dx = 1. We then get an 
averaging operator: 

Let Pu;,T be the associated probabihty measure: 

P.,t(/ eA) = ^J^ ^A{f{t))w (^) dt . 

Note that the requirement w G C^(0, oo) imphes that the Fourier 
transform of w decays rapidly: 

w{x) <^ la;]""^, as |x| — >• cxd 

for all A > 1. In the concluding section |H1 we will relax the restrictions 
on w to allow other averages, e.g. w = l[i^2] so that we take t uniformly 
distributed in [T, 2T], or w(t) = 211^ ^i when we take the eigenvalue 
A = 1/4 + t^ uniformly distributed m'[E, 2E], E = 1/4 + T^. 

6.2. The expected value of S/l- We will first show that the mean 
value {Sf^L)^j, tends to zero as T ^ oo provided L = O(logT): Aver- 
aging (j5.6|) we find 

Note that since \ogAf{n) ~ 21ogra and supp/ C [—1,1], the sum is 
over n < e'"^. Using w{x) <^ x~^ as x — > oo, we have 

Since '^n<x Pi^Y <^ a; by LemmaEHl we have by the Cauchy-Schwartz 
inequality that 

gTrL 

which goes to zero since we assume L = O(logT). 

Note that this argument also works when we allow straight averages 
(such as w = l[i,2]) as long as L < - logT. 

6.3. The variance of S'/l- 

Proposition 6.1. //limsupTrL/logT < 1 then as T —>■ oo.' 

(6.1) ((5/,.)^L,.--i=:^/°°/>)V-d. 

where k is given by ()3.7|) . 
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Note that we have 

ttL/2 

e(i-).W2 « ^^ « e 

for all e > 0. 

Proof. To compute {{Sf^i)'^)^j,, use fl5.6|) to get 

,\ogAf{m) c- logA/'(r2) 
pyiiijpyiLjj \ 

m,n<e'" 



/iQ ^2\ 1 '^ ,^^ ^/Q^ w/logAr(m) c-logAr(? 



T 



X ^ il;(— (eilogA/'(m) + e2logA/'(r2))) 

ei,e2=±l 

We now deduce that as T ^ oo, the only non- vanishing contribution 
is from the "diagonal terms" where ei = — e2 and 

J\f{m) = J\f{n) 

that is m = n. 

If m 7^ ra we may use Lemma 14.21 to get a lower bound 

(6.2) I logAr(m) ± logAf{n)\ > 



mm m,, n 



To be included in the sum, we need M{'m),M{n) < e'^'^^, that is m, n < 
e'^^, and so 

wi^ie, log ATM + e, logATH)) « ( "'^^l^-^V « C^f ■ 

Moreover, from ^„<a,/3(^)^ ^ x we get by Cauchy-Schwartz that the 
off-diagonal contribution is dominated by 

for all A > 1. This goes to zero if ttL < (1 — 6){logT) for some 6 > 0, 
which we assume. 

The diagonal terms m = n give 

^ ' n>2 

(where we used -^(0) = 1). Since there are two such terms (correspond- 
ing to ei = — e2 = +1 or —1), we have the total diagonal contribution 
being 

^ ^ n>2 
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This can be evaluated asymptotically as L — > oo using Peter's formula 
(Lemma |3.1|) to give 

Thus we find 

if limsupvrL/logT < 1. D 



7. Higher moments 

We can now show that Sf^ii^r) has Gaussian moments: 

Theorem 7.1. For K > 3 the K-th moment of Sf^i/ai converges 
to that of a normal Gaussian provided that L -^ oo with T but that 
L = o{\ogT): 



Ml 



K = 2k even 
lim 

^ ^ "''^ [O , K odd 

7.1. Reduction to the pre-diagonal. By fj5.6|) the K-th moments 
of Sf^L is given by 

(7.1) ((s„)->„,, ^ -i^ E n^(„,)/(!2^) 

ni,...,nK<e'"^ i=l 

We now show that as T -^ oo, the only (possibly) non- vanishing 
contribution to ()7.1|) is for terms satisfying: 

K 

Y,Vj^ogAf{n,) = 

j=i 
that is we have 
(7.2) 

^ ' '?j=±lE,'?ilogA^K)=Oi=l 



+ 0{ 
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for some aK,'~iK > 0. Since L = o(logT) the remainder term vanishes 
as T — > oo. 

To prove this, recall that by Lemma 1^^ if ^,=i^j logA/'(?2j) 7^ 
then for some 6k > 

K / ^ \ ~^^~^+^/^ 

i=i \j=i / 

since only terms with M{nj) < e^'^^ appear in ()7.1|) . Thus for these 
terms we have 

j 

Replacing /3(n) by log^ n ^ L^ in ()7.ip gives that the contribution of 
the terms with J2j=i Vj log A/'(nj) 7^ is dominated by 

ni,...,ni(<e^^ 

Since L = o(logT), this vanishes as T ^ 00 (in fact we need only 
assume that L < c^^logT for this, if ck is sufficiently small). This 
proves ()7.2p . 

7.2. Off-diagonal terms. In ()7.2|) we consider the sum of non-diagonal 
terms, that is terms for which there is at least one index j such that 
Uj 7^ Hi for all i ^ j- To handle these, we use Lemma f4. II which forces 
the relation 

K 

l[Ar{n,r^ = 1 

i=i 
to decompose into a union of such relations. Thus there is a decompo- 
sition 

{1,2..., K} = Y[S, 
so that in each subset Sj we have 
(7.3) Y[^f(n,r^ = l 



and the norms Af{ni) = {rii + ^Jnf — 4)/2 lie in the same real quadratic 
field Q(-^/dj) for all i E Sj. In the diagonal case there are K/2 such 
sets, e.g. Si = {1, K/2 + 1}, ^2 = {2, K/2 + 2}, . . . and the identities 
are of the form J\f{nj)J\f{nK/2+j)~^ = 1, j = 1, • • • , K/2. 
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In the off-diagonal case we assume that there is a subset Sj contains 
at least 3 elements. The number r of subsets is then at most {K — l)/2, 
since 

r 

i=i 
To count such tuples of nt, we denote for each subset Sj by dj the 
common value of the square-free kernel oinf — 4, i E Sj and then write 

n^ - 4 = djf^, I e Sj 

Let e{dj) be the fundamental unit of the field 'Q_{^fdj) and write 
Mini) = e{dj)'^^\ i G 5*^. Since logM {rii) <^ L we have ki <^ L/ \oge{dj), 
i G Sj and the relation ()7.3|) implies ^jg5, ±fcj = 0. Thus for each sub- 
set Sj there are at most 0((L/loge((ij))*'^^~^) solutions of ()7.3j) with 
\ogJ\f{ni) < L. 

Recall that we are summing over logA/'(r7,) < 27rL. Using l3{n) ^ 

(logn)^ -C L'^ we find that the off-diagonal contribution is bounded by 

the sum over all partitions {1, . . . , K} = lYj=i ^j of 

(7.4) 

'' T 
L^T\ Y {-. ttt)^*^^^'^^' < L'^md fundamental : eid) < e'^HY 

where r < {K — l)/2 is the total number of subsets Sj in our partition. 

Lemma 7.2. The number of fundamental discriminants d > for 
which t{d) <X IS 0{X^+^) for all 5 > 0. 

Proof. We need to bound the number of fundamental discriminants d 
for which the fundamental solution e{d) = {xd+ydyd)/'^ oix^ — dy"^ = 4 
is at most X. Since e{d) ~ x^, this is equivalent to bounding the 
number of fundamental (i's for which x^ <^ X. In turn, this number is 
majorized by the number v{X) of all triples (d, x, y) of positive integers, 
with d = 0, 1 mod 4, for which x^ — dy"^ = 4 and x < X, which is the 
sum 

z/(X) = '^i^{d,y> l,d = 0,l modA: dy^ = x^ -A} . 

x<X 

Since for x ^ 2 the number of pairs {d, y) with dy"^ = x^ — 4 is at most 
the number of divisors r(x^ — 4) of x^ — 4, we find that 

K^)< Y. ^(a;'-4)< Y. ^'^x^^' 

2<x<X 2<x<X 

for all (5 > 0, by virtue of the bound T{n) <^ n^ for all 5 > 0. 
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Note: A more refined argument [2S1 Lemma 4.2] shows that z^(X) 
is asymptotic to ^X, so that one can replace the bound 0{X^^^) by 
0(X). D 

Thus we find that ()7.4|) is bounded by 

for all 5 > 0. Since ai > e^^^^^"-^/^ for all e > 0, this shows that the 
sum of the off-diagonal terms is 0(0"^^"^"'''^), for all e > 0. To prove 
Theorem 17. II it thus suffices to evaluate the diagonal contributions. 

7.3. The diagonal contribution. Assume now that there is the same 
number of + signs as there are — signs. That is K = 2k is even, and 
there are ( ^ ) such choices of signs. For simplicity assume the first k 
are + and the last k are — . Thus we have to evaluate the sum 

n-=iA/'K)=n|'fc+iA^K)^=^ 

There are k\ ways to pair off variables from the first k and the last k, 
such as the pairing rij = Uk+j, I < j < k. Each such pairing contributes 
a term 

There are overlaps between the different ways of pairing off variables, 
which correspond to intersection of diagonals such as rti = n2 = n^ = 
71,4. The contribution of these was already estimated in the study of the 
non-diagonal terms, as they correspond to relations ()7.3|) where some 
subset has more than two elements. 

Thus the total contribution of diagonal terms to ((5'/^l)^'^)^ j, is 
asymptotically 

2fc\ M ,^i,fe (2A;)! 2fc 
This proves Theorem 17.11 D 



8. Conclusion 

Since the Gaussian distribution is determined by its moments. The- 
orem 17.11 implies 
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Theorem 8.1. Assume that L —>• oo as T —>• oo but L = o(logT). 
Then 



lim P^T ^^^^^^^ <x)= e~" /^ 



'r-"^ ' ^L J-oo V27r 

So far we have assume that the weight function w defining the aver- 
ages is in C^(0, oo). To deduce the resuhs for the standard averages 
{w = l[i,2]) as in Theorem 12. 11 one proceeds by approximating l[i 2] by 
"admissible" w's in a standard fashion, see e.g. JH]- The details are as 
follows: Fix e > 0, and approximate the indicator function l[i,2] above 
and below by smooth functions x± > so that X- < l[i,2] < X+) where 
both x± aiid their Fourier transforms are smooth and of rapid decay, 
and so that their total masses are within e of unity: | Jx-|-(a;)da:; — 1| < e. 
Now set u!± := X±/ f X±- Then u!± are "admissible" and for all t, 

(8.1) (l-eV_(t)<l[i,2](t)<(l + eV+(t) 

Now 

SfAr) ^ A_ r. fSf,L(T)\^ (t 



meas jt G [T, 2T] : ^^ ^ ^\ ^ J ^"^ l^T^j ^''''l [t > "^^ 
and since ()8.ip holds, we find 



;i -e)Pa._ T i ■^' ' ' G ^ ^ < - meas a G [T, 2r] : -^'^^ ^ G ^ 

^/,L(r) 



< (1 + e)P^^,r ^^^^^ G ^ 

By Theorem 18. II we find that 

(1-e)^ [ e-^'/2 dx < lim inf- meas It G \T,2T] : ^^Miill ^ ^j 
V2^y^ - T^oo T I (Tl J 

with a similar statement for limsup; since e > is arbitrary this shows 
that the limit exists and equals 

lim 4 meas (t G [T, 2T] : ^M^ll e a\ = -^ [ e^^'/^ ^^ 

which proves Theorem 12. II 

The same consideration applies to other positive statistics, such as 
the number variance. 
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